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Abstract 

We point out new ways to search for charmless baryonic B decays: baryon 

pair production in association with rj' is very likely as large as or even a bit 

larger than two body K-k/tttt modes. We extend our argument, in weaker 

form, to B ^ 'J + Xg and ii' + X. Although calculations are not reliable, 

estimates give branching ratios of order 10~^-10~^, where confidence is gained 

from recent experimental finding that B D*pn, D*ppTT are not far below 

D*Tr and D*p rates. Observation of charmless baryon modes would help clarify 

the dynamics of weak decays to baryonic final states, while the self-analyzing 

prowess of the A-baryon can be helpful in CP- and T- violation studies. 
PACS numbers: 13.25Hw, 13.20He, 13.60Rj, 14.40Nd 
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Many charmless mesonic B M(^s^M decays have emerged at the 10~^ level in recent 
years, giving evidence for strong b —>■ s penguins and tree level b —>■ u transitions. In 
contrast, the search for rare baryonic modes has been less fruitful. The most recent limits 

are B Ap, Apvr" and pp < 0.26, 1.3 and 0.7 xlO~^, respectively, improving previous 
bounds 0] by more than an order of magnitude, except for the pp mode which has a 2.8a 
excess. Theoretical work on rare baryonic (we use M and B to stand for meson and baryon, 
respectively) decay is equally sparse 1^-^, but in general they predict B B(s)B to be below 
10^^, oftentimes considerably below. With the advent of B factories, two body baryonic 
modes should eventually emerge. But one may wonder: Where is the best place to search 
for charmless baryonic modes? 

We suggest that charmless baryon-antibaryon final states in B decays may show up in 
association with t]' and/or 7 with sizeable branching ratios, i.e. ~ 10^^ — 10^^ or more. 
Although theoretical calculations are very unreliable, the bright side of this is that theory 
will learn much from experiment once the measurements become available. In particular, we 
would get important input for understanding the dynamics of weak decays. Furthermore, 
baryonic final states offer new observables that should be sensitive probes of CP and T 
violation. 

We take cue from the surprise discovery of the B ^ rj' + X modes. Without guidance 
from theory, the CLEO Collaboration discovered that both inclusive B —>■ rj' + Xs (where 
Xg = K + mT) and exclusive 5 ^ rj'K modes are very large (> 6 x 10~^ and ~ 8 x 10~^, 
respectively). To date they are still not clearly understood. Theoretical work done after 
the fact 1^ can still not fully explain the exclusive rate even by some ad hoc tuning of 
parameters. For the inclusive mode, an interesting but still controversial proposal |[TO| , pT| , 
based on 6 sg* followed by g* gr]' transition (which is motivated by the gluon anomaly), 
seems to account for the observed mx^ recoil spectrum and rate. In the following, we 
give two semi- quantitative arguments, one from the inclusive perspective with the anomaly 
mechanism, the other from exclusive perspective and based on a pole model, that suggest 
B — > ?7'BsB(7r) |]l2l may be comparable to the rj'K modes with Br ~ 10"^. Interestingly, 



our suggestion is supported by the very recent observation [|T3| of S ^ D*pn and D*ppTi 
modes with rates not far below B D*7r, D* p. 

The partial reconstruction technique used in inclusive B rj' + Xs study was initially 
developed for inclusive b ^ S'y studies. Since both the inclusive 6 — > 37 ~ 3 x 10~^ 
and exclusive B K*'~f ~ 4 x 10^^ |TB] are comparable to their rj' counterparts, we extend 
the pole model argument and suggest that B 7BsB(7r) may also be promising. Similar 
arguments suggest that B iuBB^n) should also be searched for. We also comment on 
the special case of 5 — > J/ipBsB |T^, where phase space is extremely limited. 

Let us understand why B — > B(s)B(7r) modes are so suppressed. Baryon formation in B 
decays is more difficult than the mesonic case, as reflected in model calculations. In pole 
models 0,0, the strong B B^Bi transition (B;, is a 6-baryon) is followed by a B;, B2 
weak transition, where the large imbalance in mass is the main source of uncertainty. The 
estimate of transitions to final states involving spin 3/2 baryons |^ seem to be ruled out 
by experiment already Another intuitive picture involves diquarks, which we denote 
generically as D. One first has a 6 — Dig weak decay |^, which, together with the spectator 
quark, gives a D1D2 pair. Further creation of a qq pairs leads to a BB final state. However, 
since the energy release in B decays is so much larger than the argued color 3 diquark 
binding scale ~ 1 GeV, the approach is dubious for charmless two body final states. Finally, 
the QCD sum rule approach tries to evaluate directly the i?-B(s)-B three point function, 
and it is the only method that has studied penguin effects so far. The basic problem may 
be the applicability of sum rules to B decay to light hadrons. Analogous to the diquark 
picture, one relies on a soft q-q pair creation model, hence it seems to be better suited for 
B — > BcB processes where the energy release is lower. 

We can now understand why B — > B(s)B modes are suppressed compared to B ^ M(^s^M: 
Baryons are more complicated objects than the "atomic" mesons and hence harder to form. 
The weak Hamiltonian induces b — > dqq and sqq transitions that lead to final states already 
containing two qq pairs, typically with matching color that project easily onto M(s)M final 
states. In contrast, not all the ingredients are present for the B — >■ B(s)B case. The need for 
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creating another qq pair leads to suppression by either a strong couphng, or by the intrinsic 
softness of nonperturbative pair creation against the rather hard weak decay. 

Applying the diquark model to penguin processes serves to illustrate further the point. 
The penguin process can be viewed as induced by a 6 — * sg* transition where g* is virtual. 
One may think that gf* — DD can be treated on the same footing as g* — qq, hence 
6 — s> sDD could be analogous to 6 ^ sqq. However, the picture fails because, while quarks 
are truly fundamental, diquarks are quark-quark correlations at best up to typical hadronic 
scale. Since the g* virtuality in two body penguin transitions is well above this scale, the 
g* — i> DD transition is suppressed by some form factor. Thus, the smallness of -B — > B(s)B 
modes is rooted in the large energy release. They are further suppressed compared to 
B M(s-)M modes because of the more complicated composition of baryons. 

From discussions above, we have gained some insight into where charmless baryonic B 
decays may be larger: One has to reduce the energy release and at the same time allow for 
baryonic ingredients to be present in the final state. A natural starting point is the inclusive 
B ^ rj' + Xg decay, where a large rate of ^ 6 x 10"^ is observed for p^ji > 2.0 GeV. Much 
energy is already carried away by the rj' while the signal is established by requiring a cut 
1^ on recoil system mass nix^ < 2.35 GeV. The observed inclusive mx^ spectrum is so 
far accounted for only by the anomaly mechanism [|l^,|ll|], namely b sg* followed by 
g* —>■ grj' with effective coupling motivated by the gluon anomaly. It has been argued that 
the anomaly coupling should be form factor suppressed since the g* is rather virtual 



[V(p ~ 3 GeV). However, the problem is interestingly nontrivial |jTT| because of high gluon 



binding scale in the Gf^yG^^^ channel, which has no analog in the 7* 77r case. At any 
rate, we take this as a model that is rather effective in producing fast rj' in B decays. We 
therefore consider the transition bq —>■ r]' + sgq where sgq forms a color singlet. 

Treating the gluon as a parton in the final state, the sgq system gives an mx^ recoil mass 
spectrum that is in good agreement with data ^ and peaks roughly at 2.3 GeV. Although 
sgq clearly can evolve into K + nir, it is instructive to visualize how it may feed a single 



kaon. In Ref. |Tl[ an effective nig = 0.5 GeV was used in final state phase space to remove 



"soft" gluons (below consituent qq threshold). Since there are no infrared singularities, it 
was pointed out that the mx^ region covered by nig ^ 0.5 GeV might be swept under the 
kaon, and could by itself account for the observed size oi B ^ rj'K. Such ^^\sgq) Fock 
component of the K meson" contributions to -B ^ 1]'K have not been taken into account in 
the traditional approach [0. Here we extend this observation and exploit ^f* DD diquark 
pair creation to construct baryonic final states. 

We illustrate the bq ^ r]' + sgq transition and sgq sDDq evolution in Fig. 1(a) and 
(b), respectively. Although the g-D-D coupling can be quantified, we shall refrain from 
introducing further model dependence in baryon formation from, say, Dg. Rather, we shall 
utilize mainly kinematic arguments. Following the example of how a single K may emerge 
from sgq, we take nig as an effective mass in final state phase space to correspond to qq or 
DD formation thresholds. For nig ^ 0.6 GeV ^ 2nig, where nig is the constituent quark 
mass, one can view that the gluon "has no place to go" , hence it could only end up in the 
K meson, as argued earlier. For 0.6 GeV ^ nig ^1.1 GeV ^ 2nis ~ 2mQ (we treat and 
mQ as roughly equal), the "gluon" can only split into uu and dd, and one has non-resonant 
formation of Kn, K2'k, etc., or the formation of \Kg) = \sgq) hybrid mesons. For ^ 1.1 
GeV, not only uu, dd and ss are open, so is DD. Until effective becomes very massive, 
say beyond 1.8 GeV, diquark pair formation is on similar footing with qq and is not form 
factor suppressed. 

We depict in Fig. 2 the regions separated by nig = 0.6, 1.1 and 1.8 GeV. Counting spin 
degrees of freedom only, we estimate that roughly 1/13 and up to to 3/7 (depending on scalar 
or vector diquark nature) of the rate in the 1.1 GeV ^ nig ^1.8 GeV domain corresponds 
to sDDg final state. For nix^ < 2.3 GeV (the recoil peak in anomaly model), phase space 
and kinematics considerations suggest that this final state would preferentially end up in 
two body BgB final states, given that the Ap threshold is at 2.05 GeV. It is likely that one 
would receive threshold enhancement since diquark pairs are already produced to the right 
of the nig ~ 1.1 GeV curve. The modes to search for are therefore B rj'AN and similar 
low lying B^B states accompanying a relatively fast rj'. We stress that the reconstruction of 
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rj'BsB modes should be easy and with httle background, since the A+iV threshold does not 
open up until 3.22 GeV. It therefore may offer an important probe into the higher mass mx^ 
spectrum not afforded hj Xs = K + nn modes. One can in principle pursue the inclusive 
study of -B ^ 7]' AN + mr {'-/), where 7 could come from e.g. S — > A7. 

The picture we outline above bears some similarity with the explanation for the low pj/^ 
"bump" in the feed-down subtracted inclusive primary J/ip momentum spectrum seen by 
CLEO |]18|. The excess for pj/^ < 0.6 GeV is suggested [|l6l to be S ^ J/tpAp where there 
is only 128 MeV available kinetic energy. It is known that B ^ J/ip + X decay has a large 
component coming from cc color octet. Although this excess color is not necessarily shed 
by a single gluon as in our anomaly model diagram, kinematic arguments were also used 
to argue for Ap in final state. The enhancement may come about because nonperturbative 
effects are operative for such low kinetic energy. It should be noted that, because of the 
latter, the detection of Ap system recoiling against J/ ip would not be easy. In the anomaly 
model mechanism for explaining fast rj' production in B decays as we outlined above, the Xg 
recoil system has mx^ peaked at 2.3 GeV. On the one hand this is not far above Ap threshold 
so one again does not expect the opening of many channels. On the other hand, the A and p 
baryons have considerable kinetic energy since they are recoiling against an energetic 77', the 
energy of which is (conjectured to be) fed by the g*gri' vertex. Thus, discovery oi B tj'BsB 
modes with energetic t]' may be more straightforward than detecting B J/tpAp, and it 
may also give credence to the anomaly mechanism itself. 

An alternative approach offers complementary support for our discussion above from a 
different perspective. Using naive factorization and simple pole model ideas, the B rj'Ap 
decay is seen (Fig. 3(a)) as occurring in two steps: B ri'+''K" , followed by "i^"^ Ap, 
where "i^'" denotes an off-shell kaon. The first vertex can be normalized to the observed rate 
for B Tj'K, but it is very difficult to make reliable statements about the strength of the 
dimensionless K-p-A effective coupling, Qbeff- The crude approximation of g'^^jj/An ^ 0.3 
gives r{B — > r]'Ap)/r{B — > rj'K) ^ 0.3, comparable to the estimate made above through 
diquark arguments. 
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Although the simple factorization and pole model ideas are far from reliable, fortunately, 
some very recent results from CLEO offer support for the above number. Following a 
suggestion by Dunietz [|19|, the first exclusive B decays to nucleons have just been observed, 
with [0 B° D*^pn, D*'pp'K^ ^ 1.45x10"^, 6.6x10""^, respectively. As illustrated in Fig. 
3(b), starting from B^ D*"7r+, D*^ ^ 2.8 x 10"^, 6.7 x 10^^, with effective couplings 
analogous to K-p-A case above, one easily attains order of magnitude understanding of the 
large strength of D*^NN{tt) modes. In general, CLEO's observation of sizable exclusive 
B D*NN{tt) decays that are not far below B D*tt gives strong support to our 
argument that B ri'AN{7[) may not be far lower than B —>■ rj'Kl 

The situation with regard to B ^ Apy is quite similar in the pole model picture, except K 
is replaced by K*. Similar estimates as above again give r{B Apy)/r{B K*'~f) ^ 0.1- 
0.3. These estimates place the baryonic branching ratios with rj' and 7 in the range of 
10~^-10~^ and therefore within reach of the luminosities of the 5-factories. As in the rj'Ap 
case, we stress that the final states 7Ap(7r), with energetic photon characteristic of 6 — > 57, 
are reconstructible, clean and should have little background. 

It is tempting to extend the pole model picture further into B — ^ i~^i> + BB(7r) via ir/p 
poles. In principle, similar numbers would follow. We caution, however, that, unless Ei or 
rriiu are large, the energy release may not be reduced sufficiently. Furthermore, charmed 
baryon production in semileptonic B decays, where one has reduced energy release, is rel- 
atively suppressed |20|. In particular, T{B pe+u + X)/V{B e+z/ + X) < 0.015 |20 



suggest that B i^u + BB(7r) may be less promising than in association with fast rj' or 7. 

Before we conclude, let us review, in descending order of inclusive rate, the processes to be 
studied for charmless baryonic B decays. The B J/ip + X process is at 1% level. In some 
sense it is not really a charmless final state, and has very limited phase space for X = Ap 
(the only possibility). Interestingly, a distortion or a bump in the low J/ momentum region 
is found, indicating that B —>■ J/ipAp could be of order 4% - 5% of B —>■ J/ip + X rate. 
But this mode is not easy to reconstruct because of the very slow proton. The charmless 
B — > i'^u + Xu process is at 10~^ level, but, as mentioned, a bound already exists indicating 



that charmed baryon content of S — > £'^i/ + Xc is less than 1.5% of semileptonic rate. Hence 
charmless baryon content of semileptonic decays may be less promising. The B ^ rj' + Xs 
process is at 10~^ level, and, if anomaly mechanism is correct, we expect Xg to be of order 
10% or more composed of 656(71), with supporting arguments from pole models as well as 
the newly observed B D*NN{tt) modes. Hence B r]'BsB{TT) should be searched for. 
The B J + Xg is at 3 X 10~^ level. Although one only has the pole model argument to 
suggest the possibility of sizable baryonic recoil, it again should be studied. 

It should be stressed that, besides yielding useful information for theory about the dy- 
namics of weak decays, the observation of baryonic final states suggested here could open 
up a new direction of studies. It is well known that the decay of the A self-analyzes its spin, 
which can be important for studying CP- or T- violation. In particular, use of the A spin will 
allow one to construct CP-odd, Tjv-odd observables such as s^- {PpXPA), which are driven 
by the real (or dispersive) part of the Feynman amplitude and not by the imaginary (or 
absorptive) part of the amplitude, as is the case for a CP-odd, T^r-even observable, such as 
partial rate asymmetries. Since rj' and photonic modes are dominantly loop-induced hence 
sensitive windows on new physics, such probes may turn out to be very useful. 

In conclusion, we point out that B — > r^' -1-656(71) could be the most promising charmless 
baryonic modes. The rj' should still be fast, and for mx^ that is not far above Ap threshold, 
the baryonic recoil system is simple and of low multiplicity. These modes not only could 
be the first charmless baryonic modes to be detected, their detection could strengthen the 
anomaly picture, and provide new probes for CP and T violation by bringing in the powerful 
self- analyzed A spin observable. Though the argument gets a bit less compelling, a parallel 
program should also be started to reconstruct baryonic modes in the recoil system against 
the photon in P j + Xg. The traditional search for two body B — >■ 6(^)6 modes of course 
should continue, but observation of rare charmless baryonic B decays proposed here could 
open a new program for the study of interplay of weak and strong dynamics, and offer very 
important probes of new CP and T violation observables. 
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FIG. 1. Diagrams illustrating baryon formation via (a) the bq ^ rj' + sgq transition, followed 
by (b) sgq sDDq evolution, where D denotes a diquark. 
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FIG. 2. Illustration for phase space argument of 7]'+ baryon pair formation following the mech- 
anism of Fig. 1. Solid, dash, dots and dotdash curves correspond to taking nig = 0, 0.6, 1.1, 1.8 
GeV in phase space. The "gluon mass" of 0.6 (1.1) GeV mark the opening of uu, dd {ss, DD) 
thresholds, while beyond 1.8 GeV, diquark pair formation may suffer from form factors. The dotted 
vertical line indicates the experimental cut on m = in S — ^ ry' + if + mr search, while the two 
solid vertical lines to the left and right corresponds to the AA^ and A+iV thresholds of 2.05 and 
3.22 GeV, respectively. 




B 



N. ^ n 



B 



P 



n 



P 



(a) (b) 

FIG. 3. Pole model diagrams illustrating (a) B r]'Ap and (b) B D*~np production as 
mediated by 5 — > rj'K, and B — D*~Tr, D*~p, respectively. 
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